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1. Introduction 



The explanation and description of confinement in gauge theories is known to be one of the 
most fundamental problems of modern Quantum Field Theory (see e.g. [1,2]). In general, by 
confinement one means the phenomenon of absence in the spectrum of a certain field theory 
of the physical |in) and |out) states of some particles, whose fields are however present in the 
fundamental Lagrangian. A natural conjecture here is that due to a linearly rising confining 
interaction, (anti)quarks cannot exist as free particles, but form colourless bound states of hadrons. 
The most natural quantity for the description of this phenomenon is the so-called Wilson loop 
average. In the case of Quantum Chromodynamics (QCD), this object has the following form 



which is nothing else, but an averaged amplitude of the process of creation, propagation, and 
annihilation of a quark- ant iquark pair. Here, stands for the matrix-valued vector-potential 
of the gluonic field, C is a closed contour, along which the quark-antiquark pair propagates, P 
stands for the path-ordering prescription, and the average on both sides is performed with the QCD 
action. The interaction potential of a pair of static quark and antiquark can be determined from 
the Wilson loop average, defined on the rectangular contour C, according to the Feynman-Kac 
formula, 



where T stands for the size of C in the temporal direction. In particular, for large contours C, 
the exponential dependence of the Wilson loop average on the area of the minimal surface, S min ., 
encircled by C, 



(the so-called area law behaviour of the Wilson loop), yields then the linearly rising confining 
potential 



Here a stands for the so-called string tension, and R is the relative distance between a quark and 
an antiquark. This exponential asymptotic dependence of the Wilson loop average on |S min J is 
just the essence of the Wilson criterion of confinement [3]. 

The linearly rising potential (2) admits a simple physical interpretation in terms of a string 
picture: the gluonic field between a quark and an antiquark is compressed to a tube or a string 
(the so-called QCD string), whose energy is proportional to its length, and a is the energy of such 
a string per unit length. This string plays the central role in the Wilson's picture of confinement, 
since with increasing distance R it stretches and prevents quark and antiquark from moving freely 
to macroscopic distances. In order to get an idea of numbers, notice that according to the lattice 
data [4] the distance R, at which the Wilson criterion of confinement becomes valid, is of the order 
of 1.0 fm, and the string tension is of the order of 0.2 GeV 2 (see e.g. [2]). 




(1) 



V = -lim±(W(C)), 



(W(C)) -> e 



V conL (R) = aR. 



(2) 
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As it is well known, all observable dimensional quantities in QCD such as hadron masses, and 
in particular the string tension, are proportional to the corresponding power of the QCD scale 
parameter AqcdQ- In particular, one has 



167T 2 

°" « A CD = » eX P -JTWr 2»r \ 9 7T 



(3) 



where <7qcd(a0 is the QCD coupling constant measured at a scale /i. Eq. (3) means that all the 
coefficients in the expansion of the string tension in powers of S'qcdCa*) vanish. This conclusion 
tells us that the QCD string has a pure nonperturbative origin, as well as the phenomenon of 
confinement, which leads to the process of formation of such strings in the QCD vacuum, itself. 

The problem of studying the properties of the QCD string is closely related to the problem of 
finding a string representation of gauge theories possessing a confining phase. In this respect, all 
methods of derivation of the string effective action from the action of gauge theories are of great 
importance. Some progress in the solution of this problem has recently been achieved for the case 
of QCD in Refs. [5-7] using the so-called Stochastic Vacuum Model (SVM) of the QCD vacuum 
[2,8]. 

The idea which lies behind this approach is that important quantities like Green functions can 
be expressed in terms of the Wilson loop [3,8]. The main strategy of SVM is then not to evaluate 
the Wilson loop itself, but to express it via gauge-invariant irreducible correlators (the so-called 
cumulants) of the gauge field strength tensors. Such correlators have been measured in the lattice 
experiments both at large and small distances [9], which enables one to use them in practical 
calculations of various physical quantities. In order to express the Wilson loop average (1) via 
cumulants, one makes use of the non-Abelian Stokes theorem and cumulant expansion (see [2] for 
a detailed discussion and Refs. therein), which yields 



--Itrexp _29p J da ^ x ) J da Xp (x')((F^(x)$(x,x')F Xp (x')$(x',x))) , (4) 

\ ^min. ^min. / 

where <$>(x,y) stands for the parallel transporter factor along the straight line. In the derivation 
of Eq. (4), one has used the so-called bilocal or Gaussian approximation, according to which all 
higher cumulants can be disregarded with a good accuracy due to their smallness w.r.t. the bilocal 
cumulant. This approximation is indeed confirmed by the present lattice data [9]. 

The bilocal cumulant standing in the exponent on the R.H.S. of Eq. (4) can be parametrized 
by two renormalization group-invariant coefficient functions _DQ CD and _D^ CD as follows [8] 



((F, v {x)<$>(x,x')F Xp {x'mx\x))) = J^|(5 MA ^ p - 5 mAa)jD Qcd ^ _ ^ + 

lr TIiis means that the dimensionless ratios of these quantities, e.g. the ratio of y/a to hadron masses, are 
universal (i.e. gQCD _m dependent) numbers. The aim of all the nonperturbative phenomenological approaches to 
QCD (e.g. the so-called Stochastic Vacuum Model, which will be described below) is to calculate these numbers, 
but not Aqcd itself. 



3 



+ 7^ 



2 [dx^ 



1 d 



{{x - x') x 5 up -{x- x') p 5 vX ) + ^— ((x - x') p 5 pX -{x- x') x 



v 




(5) 



Notice, that by construction the term containing the function D± on the R.H.S. of Eq. (5) 
yields a perimeter type contribution to the Wilson loop average (4). In particular, the leading per- 
turbative contribution of the one-gluon-exchange diagram to the Wilson loop average is contained 
in the function _D^ CD . 

The nonperturbative parts of the functions _D QCD and £)^ CD have been calculated on the 

lattice in Ref. [9], where it has been shown that at large distances \x — x'\, £)^ CD — x') 2 ) <C 
^qcd (( 

x — x') 2 ), and both functions decrease as e ' x X '\/ T a , where T g is the so-called correlation 
length of the vacuum, T g ~ 0.22 fm in the S77(3)-case. 

The idea realized in Ref. [5] was to consider the quantity S e s, = — lnW(£ min .) defined by 
Eqs. (4) and (5) as a starting point for the derivation of the QCD string effective action. Due 
to the smallness of the ratio (T g /R) 2 ~ 0.04, it was then possible to perform an expansion of the 
nonlocal action S e s. in powers of (T g /R) 2 , which yielded the first few local terms known in the 
standard string theory, whose coupling constants were completely determined through the bilocal 
correlator. In Ref. [6], perturbative corrections to the obtained string effective action arising due 
to perturbative gluons propagating in the nonperturbative QCD vacuum have been calculated, 
and in Ref. [7] a possible way of the solution of the problem of crumpling of the string world-sheet 
has been proposed. 

However, it still remains unclear how one can get a mechanism of integration over string world- 
sheets in QCD, since up to now the construction of the string effective action has been performed 
on the surface of the minimal area (cf. Eq. (4)). Elaboration out of such a mechanism would 
eventually completely solve the task of finding a string representation of QCD. 

The present paper is devoted to be a first step in this direction. To this end, we shall simplify 
the problem under study by considering the related Dual Abelian Higgs Model (DAHM). This 
model approach is based on the commonly accepted observation that on the phenomenological 
level, quark confinement in QCD can be explained in terms of the dual Meissner effect [10,11]. 
According to the 't Hooft-Mandelstam picture of confinement, the properties of the QCD string 
should be similar to the ones of the electric vortex, which emerges between two electrically charged 
particles immersed into the superconducting medium filled with a condensate of Cooper pairs of 
magnetic monopoles. In the case of the usual Abelian Higgs Model (AHM), which is a relativis- 
ts version of the Ginzburg-Landau theory of superconductivity, such vortices are referred to as 
Abrikosov-Nielsen-Olesen strings [12]. 

Thus, it looks reasonable to consider DAHM as a natural laboratory for probing various ap- 
proaches to the problem of the string representation of QCD (for related investigations of AHM 
see [11,13-18]). In this respect, in Ref. [14] the so-called path integral duality transformation 
proposed in Ref. [13] has been applied to the lattice version of AHM in the London limit in order 
to reformulate the partition function of this theory in terms of the string world-sheet coordinates. 
Then the same reformulation has been performed in the continuum limit in Ref. [17]. 

In this paper, we would like to demonstrate the usefulness of the path integral duality trans- 
formation for the derivation of the string representation for the generating functionals of gauge 
field strength tensor and monopole current correlators in DAHM. In particular, we find it useful to 
study a suitably extended version of this model, which includes gauge fields generated by external 
electrically charged particles called "quarks" . 
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The outline of the paper is as follows. In Section 2, we shall derive the string representation for 
the partition function of such an extended version of DAHM in the London limit. The resulting 
action will consist of two parts, the first of which will describe the short-range Yukawa interaction 
between a quark and an antiquark, whereas the second one, corresponding to the long-range 
confining potential, will contain the integral over all the surfaces encircled by the contour of the 
quark-antiquark pair. 

In Section 3, we generalize the result of Section 2 by introducing external sources of the field 
strength tensors in order to derive the string representation for the generating functional of gauge 
field correlators in the London limit of extended DAHM. In particular, we shall calculate the 
bilocal correlator and demonstrate that its long- and short distance asymptotic behaviours are in 
agreement with those found in QCD. 

In Section 4, we derive the string representation for the generating functional of monopole 
current correlators and calculate the correlator of two such currents. 

Finally, some technical details concerning the path integral duality transformation and the 
integration over the related Kalb-Ramond field are outlined in two Appendices. 

2. String Representation for the Partition Function of the Extended 
DAHM in the London Limit 

We shall start with the following expression for the partition function of the extended DAHM 



Z = J l^V^lVB^Veexpi- J d 4 x 



P \2 1 . 

- F ) + - \D $ 

/»' iiv) ~ 2 I M l 



+ a i$r-?r 



(6) 



where = e ' x ' is an effective Higgs field of "Cooper pairs" of magnetic monopoles, 

and F^y = d^B u — dyB^ are the dual gauge field and its field strength tensor, = — 2igB fl 
is the covariant derivative with g standing for the magnetic coupling constant. Notice that in Eq. 
(6), F^ v stands for the field strength tensor generated by external "quarks", defined according to 
the equation 



d F 



(7) 



with 



3u(x) = e 



dr 



dr 



S(x — x{t)) 



standing for the conserved electric current of a quark, which moves along the closed contour C, 
parametrized by the function x^ir), < r < 1, x M (0) = The electric coupling constant e is 

related to the magnetic one via Dirac's quantization condition eg = -|, where n is an integer 0. In 
what follows, we shall for concreteness restrict ourselves to the case of monopoles possessing the 
minimal charge, i.e. put n — 1. 
The solution of Eq. (7) reads 



2 Here, we have adopted the notations of Ref. [11]. 
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where S A1I/ (x) = / da ttv (x(^))S(x — x(£)) is the so-called vorticity tensor current [16] defined on the 

s 

string world-sheet E. This world-sheet is parametrized by the four-vector x M (£), where £ = £ 2 ) 
is a two-dimensional coordinate. Due to the Stokes theorem, the vorticity tensor current is related 
to the quark current according to the equation 

In particular, this equation means, that in the case, when there are no external quarks, the 
vorticity tensor current is conserved, i.e. due to the conservation of electric flux all the strings 
in this case are closed. Notice, that when external quarks are introduced into the system, some 
amount of closed strings might survive. From now on, we shall restrict ourselves to the sector of 
the theory with open strings ending at quarks and antiquarks only. 

In the London limit, A — > oo, the radial part of the monopole field becomes fixed to its v.e.v., 
|$| — > 77, and the partition function (6) takes the form 



z = J vB.ve^ ve^- exp j- J d 4 x [± (f, u - f*) 2 + £. (d,e - 2 9 b„) 



(8) 



where from now on constant normalization factors will be omitted. In Eq. (8), we have performed 
a decomposition of the phase of the magnetic Higgs field 9 = 9 smg - + 9 rcg -, where 9 smg -(x) obeys 
the equation (see e.g. [13]) 

e^x P dxd p 9 shlg ix) = 2ttE /u/ (x) (9) 

and describes a given electric string configuration, whereas 9 reg '(x) stands for a single- valued 
fluctuation around this configuration. Notice also, that as it has been shown in Ref. [17], the 
integration measure over the field 9 factorizes into the product of measures over the fields 9 smg - 
and 9 reg -. 

Performing the path integral duality transformation of Eq. (8) along the lines described in 
Ref. [13], we get 



= exp (-(27re) 2 J rf 4 xE 2 ^ J VB p J Vh^ J Vx^t) exp j J d A x 



--F 2 
4 w 



[XV ) 



(10) 



where H^ u \ = d^h u \ + d\h pv + d u h\ p is the field strength tensor of an antisymmetric tensor field 
(the so-called Kalb-Ramond field). Next, by carrying out the integration over the field B^ in 
Eq. (10), we obtain 



Z = j Vx^(0 J Vh^exp I J d l 



x 



11) 



12?7 2 ^ A 4e 2 

The details of the derivation of Eqs. (10) and (11) are outlined in the Appendix 1. 

Finally, the Gaussian integration over the field in Eq. (11) (see Appendix 2) leads to the 
following expression for the partition function (8) 
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Z = J T>Xp(£) exp < -7T 2 J da Xl/ (x) J > 
Is s 



d(T\ v {x) / da„ p (y)Dx U) up(x - y) 



In Eq. (12), the propagator of the field h pv has the following form 



where 



D 



\i/,pp' 



dS,Jx) 



*) = D { ± p {x) + D 



, (2) (x) 



T] 3 K 



ln 2 e \x 



t I S\,,5,,r, — S,,,,S 



\p,<Jup ~ VpuVXp , 



e?/ 



4tt 2 x 



2 Id 



m 



2 

2 1 



\x\ 2 



hpSvp - 5 pt ,5 X p + 



^ jKl + — (K + K 2 ) + ^K 3 
4 ar / 2\x\ 



m 
T 



'121 



(13) 



&\pXpX v -\- 5 pv x\x p 6 ^x v x p $vpXpX\ 



(14) 



Here Ki = Ki(m \x\), % — 0, 1, 2, 3, stand for the modified Bessel functions, and m = 2 is the mass 
of the dual gauge boson generated by the Higgs mechanism. Due to the Stokes theorem, the term 



da Xu {x) J dapp(y)D { ^ p (x-y) 



can be rewritten as a boundary one (see Appendix 2), which finally leads to the following repre- 
sentation for the partition function of extended DAHM in the London limit f\ 



Z = exp 



J exp 



-fdxpfdyp 

c c 



— J da, v {x) J da, v {y) 

EE 



(15) 



The first exponent on the R.H.S. of Eq. (15) leads to the short-range Yukawa potential, 

V Y uk.(R) « ^e- mR . 
a 

Notice that since quarks and antiquarks were from the very beginning considered as classical 
particles, the external contour C explicitly enters the final result. Would one consider them on 
the quantum level, Eq. (15) must be supplied by a certain prescription of the summation over the 
contours [1,11]. 

3 It is worth mentioning, that an analogous expression for the partition function has been obtained and investi- 
gated in Ref. [11]. This has been done by making use of techniques different from the ones applied in the present 
paper. 
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The integral over string world-sheets on the R.H.S. of Eq. (15) is the essence of the string 
representation of the partition function. Being carried out in the saddle-point approximation, it 
results in the last exponent on the R.H.S. of Eq. (15), where the integrals are taken over the 
surface of the minimal area, encircled by the contour C, S min (cf. Eq. (4)). This term then yields 
a rising confining quark- ant iquark potential (2). 

As it has already been discussed in the Introduction, the string tension a in Eq. (2) is nothing 
else, but the coupling constant of the Nambu-Goto term, 



where g stands for the determinant of the induced metric tensor g a b{€) = (^^(OX^tXO) °f 
the string world-sheet, a, b = 1,2. This term is the first local term in the derivative expansion of 
the full nonlocal string effective action 

i] 3 r r K-y{m\x - y\) 



\x - y\ 



^min. 

The second local term in this expansion is the so-called rigidity term [19,20] 



5 rigidity = — / ^V^VaVtOX^VtO), 



where £^(0 = -j^e ab {d a x^)){dbX u {^)) is the extrinsic curvature tensor of the string world-sheet, 

and ^ is the inverse bare coupling constant. 

Making use of the results of Ref. [5], it is possible to derive from the nonlocal string effective 
action S e g, both the string tension of the Nambu-Goto term and the inverse bare coupling constant 
of the rigidity term. The latter one occurs to be finite and reads 

1 ire 2 , , 

— = ~5-- 16 

a 8 

In particular, one can see that < 0, which reflects the stability of strings [20,15,21] 0. As far as 

ns 

cr = Tcr] 2 K (a) ~ tct] 2 In — , (17) 



the string tension is concerned, one obtains 

2 

7a 

where 7 = 1.781... is the Euler's constant, and a stands for a characteristic small dimensionless 
parameter. In the London limit, we get a ~ ^, where M = 2y/2Xr] is the monopole mass following 
from Eq. (6). This mass plays the role of the UV momentum cutoff somehow analogous to the 
inverse lattice spacing in QCD. Notice, that the logarithmic divergency of the string tension in 
the Ginzburg-Landau model and AHM is a well known result, which can be obtained directly 
from the definition of this quantity as a free energy per unit length of the string (see e.g. [22]). 
Clearly, both the string tension and the inverse bare coupling constant of the rigidity term are 
not proportional to g, which means that these quantities are essentially nonperturbative similarly 
to the QCD case (cf. Eq. (3)). 

4 Notice, that the negative sign of this coupling constant alone does not guarantee stability of strings. In the 
non-Abelian case, one should elaborate out another mechanisms to ensure this stability and get rid of crumpling 
of the string world-sheet [19,7]. 
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3. String Representation for the Generating Functional of Field 
Strength Correlators 

In this Section, we shall derive the string representation for the generating functional of field 
strength correlators in the London limit of extended DAHM. From this we shall then obtain an 
expression for the bilocal correlator and compare it with the one in QCD. Our starting expression 
for the generating functional reads as follows 



z [s a( s] = J vB^ve^-ve^- exp j- J d A x x - (f m „ - f^) 2 + ^ {d„e - 2 9 b ii ) 2 + iS^f, u j , 

(18) 

where Sp V is a source of the field strength tensor F^ = ^ep U x p Fx p , which obviously corresponds to 
the field strength of the usual gauge field Ap in the AHM. Performing the same transformations 
which led from Eq. (8) to Eq. (12), we obtain from Eq. (18) 

Z [S a p] = exp - J d 4 x {Sl„ + AiiieS p U Y.p^j 



■ J Vx^g) exp j- j d 4 xd 4 y ^7rE A „(a;) - ^Sxv(x)^ D Xv ,p P {x - y) (jr£ Mp (y) - % -Sp P {y) S j J . (19) 

Let us now derive from the general form (19) of the generating functional the expression for the 
bilocal correlator of the field strength tensors. The result reads 



Fxu{x)Fp P {y) 



1 5 2 Z[S a(3 ] 



Z[0] 5S Xl/ (x)5Spp(y) 



(S\p$v P - 5 X p5p V ) 6(x -y) + —D Xu ,pp(x - y)- 



-4tt 2 ^ ^2eLxu{x) - da a(3 (z)D a/3M (z - x)^j ^2eLp P (x) - ^ J da 7( (u)D l(t p P (u - y) 



where 



/P^(0(...)exp 



(...) 



-7t 2 / da a/3 (z) j da lC (u)D a p M (z - u) 



JVxp(£)exp 



(20) 



-7T 2 / da af3 (z) J da^(u)D af3r/( (z - u) 
s s 

is the average over the string world-sheets, and the term with the ^-function on the R.H.S. of Eq. 
(20) corresponds to the free contribution to the correlator. 

Let us next study the large distance asymptotics of Eqs. (13) and (14), i.e. consider these 
equations at |x| ^> ^. Then one can see, that due to the large distance asymptotics of the 

modified Bessel functions, the propagator Dx u ,p P (x) has the order of magnitude ^. Therefore, at 
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such distances the absolute value of the second term on the R.H.S. of Eq. (20) is much larger 
than the absolute value of the last term, provided that the following inequality holds 



v 2 m 



«1, 



(21) 



where |E| stands for the area of the surface E. In what follows, we shall restrict ourselves to the 
case of small enough 77 and g, for which this inequality is valid, and consequently the last term on 
the R.H.S. of Eq. (20) can be disregarded w.r.t. the second one. 

Following the SVM [2,8], let us parametrize the bilocal correlator of the field strength tensors 
by the two Lorentz structures similarly to Eq. (5) 



F Xu (x)F^(0) 



\d(x 2 ) + 



1 

+ 2 



d\ x ll 5 vp - x p 5 Ufl +d u [ 



(22) 



Then in the approximation (21), by virtue of Eqs. (13) and (14), we arrive at the following 
expressions for the functions D and D 1 



D M = 



m 3 Ki 

47T 2 \x\ ' 



and 



D, (x 2 ) 



m 



2tt 2 x 2 



\xl + j[ K ° + K \ 



(23) 



(24) 



In Eq. (23), we have discarded the free 5-function type contribution, since we are working at large 
distances so that x 7^ 0. The asymptotic behaviour of the coefficient functions (23) and (24) in 
the limit \x\ ^> — under study is then given by 



D 



mr 



AV2tt2 (m\x\)* ' 



(25) 



and 



D 1 



m 



-m\x\ 



2V2- 



VT2 772 \ x 2 



(26) 



For bookkeeping purposes, let us also list here the asymptotic behaviour of the functions (23) and 
(24) in the opposite case, \x\ — . It reads 



D 



m 



4tt 2 x 2 '' 



(27) 



and 



D 1 



21 1 4 ' 

7T Z \X\ 



(28) 
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One can now see that according to the lattice data [9] the asymptotic behaviours (25) and (26) 
are very similar to the large distance ones of the nonperturbative parts of the functions _DQ CD and 
_D^ CD , which parametrize the gauge-invariant bilocal correlator of gluonic field strength tensors 
in QCD (see Eq. (5)). In particular, both functions decrease exponentially, and the function D 
is much larger then the function Di due to the preexponential power-like behaviour. We also see 
that the dual gauge boson mass m corresponds to the inverse correlation length of the vacuum 
T g x . In particular, in the string limit of QCD, studied in Refs. [5-7], when T g — » while the value 

of the string tension is kept fixed, m corresponds to ^J DQC °(°l _ 

Moreover, the short distance asymptotic behaviours (27) and (28) are also in line with the 
results obtained within the SVM of QCD in the lowest order of perturbation theory. Namely, 
at such distances the function £)Q CD to the lowest order also behaves as j4 (which is just due 

to the one-gluon-exchange contribution) and is much larger than the function _DQ CD to the same 
order. Let us however stress once more, that Eqs. (27) and (28) contain only a part of the full 
information about the asymptotic behaviours of the functions D and D\ at small distances. The 
remaining information is contained in the omitted last term on the R.H.S. of Eq. (20), which at 
small distances might become important and modify the asymptotic behaviours (27) and (28). 

The above mentioned similarity in the large- and short distance asymptotic behaviours of 
the functions D and D\, which parametrize the bilocal correlator of the field strength tensors in 
DAHM and the gauge-invariant correlator in QCD, thus supports the original conjecture by 't 
Hooft and Mandelstam concerning the dual Meissner nature of confinement. 



4. String Representation for the Generating Functional of the 
Monopole Current Correlators 

In this Section, we shall present the string representation for the generating functional of the 
monopole current correlators in the London limit of extended DAHM. Such a representation can 
be derived by virtue of the same path integral duality transformation studied above. After that, 
we shall get from the obtained generating functional the correlator of two monopole currents and 
by making use of it rederive via the equations of motion the coefficient function D in the bilocal 
correlator of the field strength tensors. 

In the London limit, the generating functional of the monopole currents reads 



VB„V6 sin z-V6 reg - exp \ / d 4 x 



-(F 



T L{d,e-2gB,) 2 + J,j, 



where = —2gr] 2 (d^9 — 2gB fl ) is just the magnetic monopole current Q. 

Performing the duality transformation, we get the following string representation for Z [ J M ] 



Z [ J^] = exp 



m 

T 



d xJ Ax) 



Vx^) exp 



vr / dcr a p(z) / da lC (u)D a/3M (z - u) 



exp I 2ge X uap J d 4 xd A y 



g ( d 2 \ 
~2 £ ^\ dx a dy D ^>w( x - y > JM x ) J s(y)+ 



^Rigorously speaking, this is a current of the monopole Cooper pairs. 
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+7rE w (y) (JLD^wix - y)j Jp(x) 



(29) 



Varying now Eq. (29) twice w.r.t. J M , setting then J M equal to zero, and dividing the result by 
2[0], we arrive at the following expression for the correlator of two monopole currents 



{jp{x)jc{y)) = m 2 5[3j(x -y) + ^ex^e^ 

,( d 



d 2 



J dasc(z) J 



da xv (u) 



dx, 



-D XvM (x - 



2 dxady^ 
d 



D X u, w {x -y)+ 



-D 



MP: XV 



(y-u) 



(30) 



It is straightforward to see that the contribution of the term (14) to the R.H.S. of Eq. (30) 
vanishes, whereas the contribution of Eq. (13) to the second term in the square brackets on the 
R.H.S. of Eq. (30) can be disregarded w.r.t. its contribution to the first term, provided that the 
inequality (21) holds. Within this approximation, making use of the equation [2] 



(jp(x)j*(y)) 



d 2 



d 2 



y dx^dy^ /3a dxpdy a/ 

which follows from Eq. (22) due to equations of motion, we recover from Eqs. (30), (31), and (13) 
the expression for the function D given by Eq. (23). 

Notice in conclusion, that only the function D can be obtained from the correlator (30) due 
to the independence of the latter of the function D\. 



\D((x-y) 2 ), 



(31) 



5. Summary and Outlook 

In the present paper, we have derived the string representation for the partition function of 
DAHM extended by introducing gauge fields of external currents of electrically charged particles 
("quarks"), in the London limit. Such a representation yielded the confining and the Yukawa parts 
of the quark-antiquark interaction potential. By the local expansion of the obtained nonlocal string 
effective action, one gets, in particular, the corresponding expressions for the string tension of the 
Nambu-Goto term and the inverse bare coupling constant of the rigidity term. Those occured to 
be positive and negative respectively, which confirms the stability of strings (cf. also [20,15,21]). 
Besides that, both of these quantities are not proportional to the coupling constant, which means 
that the string nature of DAHM is of the same nonperturbative kind as the one of QCD. 

In Section 3, we have obtained the string representation for the generating functional of the 
field strength correlators. In a certain approximation (see Eq. (21)), this generating functional 
determined then the expressions of the two coefficient functions, D and D 1 (cf. Eqs. (23) and (24)), 
which parametrize the bilocal correlator. Those occured to be quite similar to the corresponding 
functions in QCD. In particular, it turned out that the large distance asymptotic behaviour of 
the obtained functions in the extended DAHM was in agreement with the existing lattice data [9] 
concerning the corresponding behaviours of the nonperturbative parts of these functions in QCD. 
We have also argued that the mass of the dual gauge boson in our approach corresponds to the 
inverse correlation length of the vacuum in QCD. These results together with SVM support the 
't Hooft-Mandelstam conjecture about the dual Meissner nature of confinement. 
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Finally, in Section 4, we have obtained the string representation for the generating functional 
of monopole current correlators in the extended DAHM. Then, by making use of the equations 
of motion, and the string representation for the correlator of two monopole currents, we have 
rederived the coefficient function D, which confirms the correctness of both approaches. 

Thus, we have demonstrated the relevance of the extended DAHM for the description of the 
string properties of confinement in QCD according to SVM and the 't Hooft-Mandelstam scenario. 
And vice versa, our results also support the validity of the bilocal approximation in SVM of QCD 
and provide us with some new insights concerning the structure of the QCD vacuum. In conclusion, 
all this shows the usefulness of the string representation of extended DAHM obtained on the basis 
of the path integral duality transformation. 

Clearly, it is now a challenge to apply this approach to the realistic non-Abelian case of QCD. 
In the spirit of this paper, it is in particular interesting to study dual QCD in the framework of 
the so-called Abelian projection method [23,18]. Work in this direction is now in progress. 
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Appendix 1. Derivation of Eqs. (10) and (11). 

In this Appendix, we shall outline some details of the derivation of Eqs. (10) and (11). Firstly, 
one can linearize the term \ (d^9 — 2gB^) in the exponent on the R.H.S. of Eq. (8) and carry 
out the integral over # reg - as follows 



j V6^ exp j-^ J d 4 x (d,e - 2gB,Y 



= j VC^V6 TCS - exp | J d 4 x 
j VC^5 («9 M C M ) exp if d A x 



-—Cl + iC^d.O-ZgB, 



^C 2 
2rf " 



iCjd^-2gB^ 



{All) 



The constraint d^C^ = can be uniquely resolved by representing in the form C M = 
\tixv\pdvh\p, where h\ p stands for an antisymmetric tensor field. Notice, that the number of 
degrees of freedom during such a replacement is conserved, since both of the fields C M and h^ v 
have three independent components. 

Then, taking into account the relation (9) between 6 smg - and E^, we get from Eq. (Al.l) 



/ 




d A x{d ti e-2 9 B lx f 



j Vh^VxJX) ex P \ J d4x 



12r] 



n 2 



(A1.2) 
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In the derivation of Eq. (Al.2), we have replaced T>9 smg ' by Dx^) (since the surface S, 
parametrized by x p (C,), is just the surface, at which the field 9 is singular) and, for simplicity, 
have discarded the Jacobian arising during such a change of the integration variable []. 

Bringing now together Eqs. (8) and (Al.2), we arrive at Eq. (10). In the literature, the above 
described sequence of transformations of integration variables is usually called "path integral 
duality transformation" . In particular, it has been applied in Ref. [13] to the model with a global 
[/(l)-symmetry. 

Let us now derive Eq. (11). To this end, we find it convenient to rewrite 



oxp \ -- j d 4 xF^j = j VG^exp^ j d^x -G\ v + iF^G ^ 



after which the ^-integration in Eq. (10) yields 



VB P exp 



d x 



= J VG^exp (- J dtx&jj) 5 (e pv \ p d p (G Xp - gh Xp - 2irieY, Xp )) 

= J VA p exp - J d 4 x {gh pv + 27rzeS^ + d p A v - d u A p Y 
In the last line of Eq. (Al.3), the constraint 



(Al.3) 



Sfiuxpd^ {G Xp - gh Xp - 2-KieL 







has been resolved by setting G\ p = gh\ p + 27rieEA P + d\k p — d p A\. 

Finally, by performing in Eq. (Al.3) the hypergauge transformation h pu — > h pv + d p \ v — d v \ p 
and fixing the gauge by choosing X p = — (see e.g. Ref. [15]), we arrive, omitting the measure 
factor X>A M , at Eq. (11). 

Appendix 2. Integration over the Kalb-Ramond Field in Eq. (11). 

Let us carry out the following integration over the Kalb-Ramond field 



Z = J Dh pi/ exp J dx\ 



-.m 



4e 2 



h„„ — i7[h.„,Yi 



pv'—'pv 



{A2.1) 



12 V 2 ^ x 

To this end, it is necessary to substitute the saddle-point value of the integral (A2.1) back into 
the integrand. The saddle-point equation in the momentum representation reads 

^5 {p 2 hT(p)+PxP,h^ip)+p pPu h^ip)) + ±htf'(p) = -inE vX (p). 
This equation can be most easily solved by rewriting it in the following way 



(p 2 Px», a p + m 2 lxu,«p) h^-(p) = -2m V 2 lZ Xu (p), 
where we have introduced the following projection operators 



(A2.2) 



3 For the case when the surface E has a spherical topology, this Jacobian has been calculated in Ref. [17]. 
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and 



with V u 



These projection operators obey the following relations 



\-fiu,\p 



-1 



vp,Xp 



-1 



^-pu,Xp^-Xp,af3 

(the same relations hold for P^Xp), and 

Ppu,Xp (1 - P, 



pu,pX *-Xp,pv, 



'-pv,a(3 



Xp,aj3 



0. 



(42.3) 
(42.4) 

(42.5) 



By virtue of properties (A2.3)-(A2.5), the solution of Eq. (A2.2) reads 



h e xT(p) = 



2nirf 



p2 _|_ m 2 



1 + 



P 



WT 



P) 



which, once being substituted back into partition function (A2.1), yields for it the following 
expression 



cxp- 



-'V/ 



d A p 



1 



(2vr)V + m 2 



p 



m 



S^(-p)S a/3 (p) 



(42.6) 



pu,afS 



Rewriting Eq. (A2.6) in the coordinate representation we arrive at Eq. (12). 

Let us now prove that the term proportional to the projection operator 1 — P on the R.H.S. of 
Eq. (A2.6) indeed yields in the coordinate representation the boundary term, i.e. Eq. (14). One 
has 



p 2 {l - P)xu,a/3 = -ji^upPxPa + $XaPvP(3 ~ ^votPxPfi ~ hfiPuPa) 



(42.7) 



By making use of Eq. (A2.7), the term 



under study, after carrying out the integration over p, reads 

Acting in Eq. (A2.8) straightforwardly with the derivatives, we arrive at Eq. (14). However, one 
can perform the partial integration, which gives the boundary term 



(42.8) 



e 2 m 



dXf, <b dy^ 



K\{m \x — y\) 
\x - y\ 



c c 

in which one can recognize the argument of the first exponent standing on the R.H.S. of Eq. (15). 
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